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NOMENCL STORE 



b 

c ? 

a 

dg 

f, (x) 
f 2 (x) 

f 3 (X) 

g 

Gr 

h 

k 

H 

kx 

m 

Nu 

Pr 

g 

Ra 

Re 

t , 



height of spine (m) 

specific heat of surrounding fluid (J/kg*K) 

diameter of spine (m) 

differential amount cf heat (w) 

cross sectional area function for spine (m ) 

profile function for spine (m) 

perimeter function fcr spine (m) 

gravitational acceleration (m/s 2 -) 

Grashof number (dimensionless) 

heat transfer coefficient (w/m a - K) 

thermal conductivity (w/m*K) 

thermal cond. of surrounding fluid (w/m*K) 

thermal cond. of insulation (w/m* K) 

thermal cond. of spine material (w/m-K) 

fin heat transfer parameter (1/m) 

Nusselt number (dimensionless) 

Prandtl number (dimensionless) 
heat flow into spine at base (w) 

Rayleigh number ( dim ensionless) 

Reynolds number (dimensionless) 
temperature of spine (fin) surface (K) 
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^ S f * C© 


temperature of surrounding fluid (K) 


Z(tf) 


enhancement factor a function of the roll angle 


GREEK 


SYMBOLS 


Q> 


a combination of variables defined by egn (A. 9) 


P 


volumetric coefficiert of thermal expansion (1/K) 


S 


diameter of spine (m) 


t 


roll angle (degrees) 


e 


temperature excess ( K) 


% 


temperature excess at spine base ( K) 


©a 

1 


temperature excess at spine tip ( K) 

dynamic viscosity of surrounding fluid (kg/m.s) 

kinematic viscosity cf surrounding fluid (m^/s) 
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I. introduction 

Electronically steered phased array antennae frequently 
exhibit degraded performance due to the lack of proper heat 
transfer between the interior crystalline materials and the 
environment. This causes higher than optimum temperature 
environments for the crystalline materials and a subsequent 
loss of reliability[ 1 ]. 

One possible solution to such a problem considers a 
surface augmentation through the use of pin fins (or spines) 
to assist in the dissipation of the heat. 

The problem is more severe when the installation is 
stationary so that the dissipation from the spines is in the 
natural convection mode and, because of this, the focus of 
this investigation is in the area of natural convection heat 
transfer and methods of enhancenent in this domain. 

Convective heat transfer from horizontal cylinders has 
been the subject of extensive studies which were summarized 
by Morgan[ 2 ]. In addition, studies of heat exchanger tube 
bundles are summarized by Incropera and and DeWitt[3]. 
However, after a considerable search of the literature, only 
a small amount of experimental data and analytical correla- 
tions dealing with natural convection in an array of hori- 
zontal cylinders, has been found. 

Intuitively, one feels that the orientation of the array 
must have an effect on the heat transfer but apparently no 
known experimental correlations have been formulated. Thus, 
the focus of this thesis is in the effect of roll angle on 
the heat transfer rates from an array of cylindrical spines. 

Additional impetus for investigations in similar areas 
with cylindrical spine arrays derives from the electronics 
industry where such cooling arrangements are used in the 
external cooling of dense integrated circuit arrangements. 
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Dimensionless correlations using dimensionless parame- 
ters are usually developed for presentation of experimental 
results and this allows for possible application to other 
systems. In developing the solution for the heat transfer 
across a boundary layer, it has been shown by 
Schlichting[4 ], and others [3,5], that the Nusselt number 
(Nu) is a proper dimensionless parameter for natural or 
forced convection. In forced convection, the Nusselt number 
is a function of the Reynolds number (Re) and the Prandtl 
number (Pr) . However, in pure natural convection where the 
buoyancy is the only driving force, the fluid velocity is 
determined entirely by the guantities contained in the 
Grashof number (Gr) which is used instead of the Reynolds 
number to establish dynamic similarity. The Grashof number 
is developed quite thoroughly by Kreith and others [3,5] and 
quickly appears in a dimensional analysis. 

In addition, for gases the Prandtl number (Pr) is nearly 
constant so that the Nusselt Modulus (Nu^) can be considered 
to be a function of the product of the Grashof number and 
the Prandtl number. This product is defined as the Rayleigh 
or number (Ra) and is discussed to further detail in the 
following paragraphs. 

In order to establish the heat transfer coefficient and 
to develop a correlation of the data with dimensionless 
parameters, relationships for the heat transfer coefficient 
and the Nusselt number were established and measurements 
were taken which permitted its calculation. A summary of 
the derivation is provided in Figure 1.1 and a detailed 
derivation is provided in Appendix A. The relation holds 
for the condition that no heat transfer occurs at the fin 
tip. Such a condition can be effectively achieved by insu- 
lating the tip. 

After calculating the heat transfer coefficient the 
Nusselt number was then determined from the relation: 



e(x) 

at x 
at x 




h = 



k d 



s 



4 



arccosh 



b 




2 



Figure 1.1 A Relation for the Heat Transfer Coefficient 




( 1 . 1 ) 



where Nu^ is the Husselt number based on the spine diameter. 

The Rayleigh number was determined from the product of 
the Grashof number and the Prandtl number which were calcu- 
lated from the following: 



Gr 



g B(T - T ) d* 
s 



C P 

Pr = 



12 



P.a = GrPr 



( 1 . 2 ) 



where again the subscript "d" infers a dependence on the 
spine diameter and where the symbols are as defined in the 
nomenclature list. 



II. EXPE RIMENT AL APPARATUS 



Figure 2. 1 displays the configuration of the experi- 
mental apparatus and the techniques used in data acquisi- 
tion. The following paragraphs provide a more indepth 
description of the component parts. 

A. TEST DISK 

The disk, shown in Figures 2.2 and 2.3, was fabricated 
entirely of 2024-T4 aluminum plate, 6.35 cm in thickness. A 
circular disk of diameter equal to 22.86 cm was cut from 
this plate and fitted with an array of 16 cylindrical spines 
as the heat transfer surfaces. Six equally spaced radial 
holes were drilled around the perimeter of the disk for 
insertion of six 200 watt cartridge heaters which provided 
the heat input necessary to elevate the temperature of the 
spines. The heaters were powered by a controllable power 
source in order to attain a steady-state condition for data 
measurements. The disk was fully insulated with soft rubber 
refrigeration insulation (kj = 0.15 w/m-K) to minimize heat 
losses from the disk surface due to both radiation and 
convection. 

B. CYLINDRICAL SPINES 

The spines, shown in Figures 2.3 and 2.4 were similarly 
constructed of the 2024-T4 aluninum. Two sizes of spines 
were used in data gathering (0.9525 cm diameter x 3.75 cm 
length and 1.27 cm diameter x 5.08 cm length) to provide a 
wider range of Rayleigh Numbers. Extension pieces, also 
made of the same aluminum, were necessary to extend the 
actual base of the heat transfer surface (the cylindrical 
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Figure 2.1 Sketch of the Test Apparatus and Configuration 
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Figure 2.2 Test Disk and Spine Assembly 



spine) to the outer level (exterior surface) 

insulation . 



of the 



C. THERMOCOUPLES 

All thermocouples were fabricated from copper-constantin 
(Type-T) of AWG No. 30 wire. They were butt welded with a 
spherical bead of approximately 0.05 cm diameter and 
smaller. All interconnecting cables were also of copper- 
constantin in order to minimize the uncertainty and to 

preclude, as much as possible, spurious thermoelectric 

effects. All thermocouples were connected to the ACUR2X 

Autodata Nine data acquisition system for temperature 

record i ng. 



1 6 




Figure 2.3 Sketch of Front View of Test 



D. ACOREX ADTODATA NINE 

The thermocouple modules 
system’s thermocouple modules i 
tion network and provide a 
printout for up to 100 chann 
internal ice-point ccmpensatio 
baths or reference junctions we 
bration checks were made at v 
expected temperature range in 
correction factors for each t 
cabling and the Autodata Nine c 
tion, the complete disk and sp 
ice-bath in order to establish 



Disk 



in the data acquisition 
nclude an ice-point compensa- 
utomatic data scanning and 
els of input. Due to the 
n networks, no external ice 
re required. In-place cali- 
arious temperatures over the 
order to develop temperature 
bermocouple, its connecting 
ircuitry as a unit. in addi- 
ine array was immersed in an 
one of the check points. 
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Figure 2.4 lest Disk and Spine Arrangement 



£. ENVIRONMENTAL CHAMBER 



The test system 
in which the ambient 
range between 35 and 



was housed in an environmental chamber 
temperatures could be maintained over a 
135 degrees Fahrenheit. 



1 8 



III. EXPERIMENTAL PROCEDURE 



A. DATA GATHERING 

The experimental procedure was to set the disk at the 
desired roll angle using a liquid bubble level indicator and 
protractor to measure the angle. The temperature control 
thermostat of the environmental chamber was then set to the 
desired temperature. Steady state environmental conditions 
were attained in approximately six to eight hours. Next the 
heaters imbedded in the disk were energized to elevate the 
temperature of the disk and base of the spine to a tempera- 
ture chosen for a predetermined Rayleigh number. The input 
voltage for the heaters was varied and the temperatures 
monitored until a steady state condition was reached. This 
was indicated by constant input power and constant tempera- 
tures. The time to attain a total steady-state condition 
was approximately three to five hours in addition to the 
time required for the chamber to stabilize. 

A steady state condition was assumed to exist when no 
change in temperature occurred for four consecutive tempera- 
ture readings taken at fifteen minute intervals. These 
temperatures were then used for the input data for the 
calculations. Typical dataruns are presented in Appendix B. 

After completion of data taking at a set condition, 
changes were made in the roll angle, the temperature of the 
system or temperature of the environment, and the process of 
attaining and maintaining a new condition of steady state 
was again started. 
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B. DATA DEDUCTION 



Temperature data printed by the Autodata Nine during the 
individual runs was used as input parameters (corrected from 
the calibration curve for each thermocouple) for calculation 
of the desired parameters for plotting and evaluation. 

All desired output parameters were calculated by the 
computer using appropriate relationships for the thermophy- 
sical properties [3] as functions of temperature and meas- 
ured temperatures as input paraneters. The computer program 
used for the calculations is provided in Appendix E. 

Graphical display of the data was also utilized to aid 
in data analyses. These curves are shown as Figures B. 1 
through B.7. 
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IV. RESU LTS AND CORRELATION 



A. RESULTS 

During each data collection run the temperatures at the 
base and tip of each spine, as well as four environmental 
temperatures were recorded. These temperatures were used as 
inputs for calculation of the various output parameters for 
each spine. These sixteen calculations were then averaged 
to provide the data point for that set of conditions. Forty 
seven data points were obtained. Typical data runs are 
included as Appendix C. 

At a specified angle, the natural logarithm of the 
Nusselt number was plotted as a function of the natural 
logarithm of the Rayleigh number for each of the data points 
at a specified angle. These plots are included as Figures 

B. 1 through B.5 in Appendix B. This graphical display was 
then used in order to facilitate correlation of the test 
data . 

The plotted results indicated that a small dependence on 
roll angle did exist. Figures B.6 and B.7 were plotted to 
show these effects more specifically. Figure B.7 was 
plotted only for Rayleigh numbers between 579 and 586 
because of the narrowness of the overall Rayleigh number 
range obtained. 

B. CORRELATION 

For each of the curves showing log (Nu) as a function of 
log (Ra) , a least sguares curve fit routine was used and a 
linear relation established to approximate these curves 
following a method presented by Kraus and Bar-Cohen[ 1 ]. A 
correlating equation was obtained for each of the roll 
angles (0), where Z { 0 ) defines a function of the roll angle. 
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The values or 2 (^) are tabulated ia Table 



1 . 



TABLE 1 



Enhancement 


Factor 


Roll 


Enhancement 


Angle ( 0 ) 


Factor ( Z( 0) ) 


0° 


0.5916 


10° 


0.5986 


20° 


0.6075 


30° 


0 .5990 


45° 


0 .5928 



The values of Nusselt numbers obtained from the experi- 
mental results are a function of a giver. Payleigh number and 
roll angle and may be correlated via: 



0.233 

Nu. - Z( 0) Ra 
d 



(4.1) 



However, observe that the dependence on the roll angle is 
weak (Table 1) and no mathematical relationship was deter- 
mined for the values of 2 {/) . 

Results were also compared with several previously 
published cor relations for the s ingle horizontal cylinder. 
Churchill and Chu[6] suggested an eguation of the fora 



Nu d = ^0.60 + 



0.387 Ra 



( 1 / 6 ) 



• 


f0.S59\ (9/16) l 


1 + 


[p r ) j 



(8/27) 
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which correlated experimental results of several researchers 
very well. 

The values obtained in this experimental program agreed 
favorably with these predictions with an average of 19% 
deviation throughtout the range investigated. Noting that 
the effect of heat transfer by radiation were also present, 
the deviation is attributed to radiation losses primarily. 
These losses would appear to net have a detrimental effect 
on the observed roll angle dependence. 
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V. CONCLUSION AND EECOMMENDATIONS 



A. CONCLUSION 

The results of the experimental program indicated that a 
dependence on roll angle did exist. This dependence is 
relatively minor (i.e., on the order of 3 %), however for 
temperatures in the range of 200 degrees Fahrenheit this 
could mean a difference of 5 to 6 degrees in the surface 
temperature of an operating system. Therefore, further 
investigation at higher Rayleigh numbers may be deemed 
appropriate. 

B. RECOMMENDATIONS 

The following recommendations are provided for possible 
follow on projects of a similar nature. 

1. Conduct the experimental pregram with the goal of deter- 
mining and removing the radiation heat transfer. 

2. Investigate other devices and methods of temperature 
measurement in order to remove more of the uncertainty 
associated with the use of thermocouples. 

3. Conduct a similar experiment for the pitch angle. 

4. Attempt a similar program with water as the convecting 
medium. 
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APPEND II A 

DEFINITION OF THE HEAT TRANSFER COEFFICIENT 



Fins and spines (pin fins) of various geometries and 
thermal conductivities respond differently to identical and 
uniform heat sources and sinks- Similarly, there are 
numerous ways in which the temperatures and heat transfer 
coefficients of sources and sinks may vary. Important to 
the analysis of fin geometries are the constraints or 
assumptions which are employed to define and limit the 
problem and often to simplify its solution. The analysis of 
fins and spines employ the assumptions of Murray[7] and 
Gardner[8]. These limiting assumptions are: 

1. The heat flow in the fin and its temperatures remain 
constant with time. 

2. The fin material is homogeneous, its thermal conduc- 
tivity is the same in all directions and it remains 
constant. 

3. The heat transfer coefficient between the fin and 
the surrounding medium is constant and uniform over 
the entire surface of the fin (spine) . 

4. The temperature of the medium surrounding the fin is 
uniform. 

5. The fin thickness is so snail compared with its height 
that temperature gradients across the fin thickness 
may be neglected. 

6. The temperature at the base of the fin is uniform. 

7. There is no contact resistance where the base of the 
fin joins the prime surface. 
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3. There are no heat sources within the fin itself. 

9. The heat transferred through the outermost edge of 
the fin is negligible compared with that leaving the 
fin through its lateral surface. 



10. Heat transfer to or from the fin is proportional to 
the temperature excess between the fin and the sur- 
rounding medium. 

In addition to proposing a profile function for longitu- 
dinal fins Gardner[8] proposed a profile function for 
spines: 



f 2 (x) = 




(A . 1 ) 



With the appropriate value of 
for the development of a gener 
for spines. 

Figure A. 1 shows a spine o 



be seen that the spine cross : 


heat, the confining 


profile, ; 


are perfectly arbitrary functii 


tip of the spine. 


A differeni 


ture excess, 0 = t 


- t , may 


heat conduction into 


and out < 


cross section f (x) : 




dq 


= k 33? [ 



n, equation A.1 may be used 
alized differential equation 

f arbitrary profile. It can 
ection normal to the flow of 
nd the perimeter of the spine 
rs of the distance x from the 
ial equation for the tempera- 
te written by considering the 
f the element dx through the 



^(x) 



d9 

dx 



dx 
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Figure A. 1 Spine of Arbitrary Profile 

The difference in heat flow into and out of the element 
must equal that dissipated by the surface of the spine. If 
dissipation occurs 1 y convection and h is the convection 
coefficient 

dq = h f (x) 0 dx 



where f ? (x) defines a perimeter function P(:<) which depends 
upon the distance x from th < -> origin of the coordinate 
system. When convection and conduction are equated. 






hf x)6 



Upon rearrangement the generalized differential equation 
becomes 



f ^ ( x) 




df ? ( x) dO 
dx dx 



h 

— f,(x) = 0 
k 6 



(A. 2) 
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The relationship between f, (x) and f 2 (x) is 
f^(x) = 7T p 2 (x)J 2 



where f 2 (x) was defined by equation A.1. This equation A. 2 
may be written 



[f 2 (x )] 2 d2 « + 


f 2 ( x )' 


A w ^ J . . 





‘dO 

dx 



2h f 2 ( x ) 9 



= 0 



(A. 3) 



where 



f 3 ( x> = 27 Tf 2 ( X ) 



Equation A. 3 is a second-order differential equation 
with variable ccetf icionts except when the spine cross 
section normal to heat flew is constant. It may be solved 
by termwise comparison with the generalized Bessel equation. 

For the cylindrical seine, shown in Figure A. 2 



and 



d_ 

dx 




0 



For this case the exponent in equation A.1 has a value of 
zero cor res pond in ] to n = 1/2. Actually £ in equation A.1 

is replaced by d, the spine diameter. Using this value of f 2 
the general differential equation for the temperature 
excess becomes 
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Figure A. 2 Cylindrical Spine 



d 2 e 




2 

m e r 



0 



(A. 4) 



i 

where m = (4h/kd)* . 

The general solution to equation A. 4 is 



9(x) = C, e mx + c e" mx (A * 5) 

1 ^2 

where the arbitrary constants are determined from the 
boundary conditions 



9(x = b) = (A. 6a) 

9( x = 0 ) = 9 

a (A. 6b) 
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This then leads directly to the equation for the temperature 
excess profile 



9 ( x ) = 9 ££^L_mx 

b cosh mb 



(A. 7) 



and the heat flow entering the fin (spine) base is 



% 



■£— kd 2 m0, tanh(mb) 
4 b 



(A. 3) 



The value of h may be determined from equation (A. 7) by 
talcing temperature readings at the tip and the base which 
are precisely "b" meters apart. 

Thus with base and tip temperature excesses, 

9 (x = b) = 6 b 9-(x=a) = 6a. 

It is easy to see from equation (A. 7) that 

e = 9b .. 

a cosh mb 



or that 



cosh( mb ) 




from which 



mb = arccosh(9, /9 ) = a 

D cl 



(A. 5) 



Then with m 



(4h/kd j't 



one may then solve for h. 



kd f a v 2 
“ 4 V 
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APPENDIX B 

GRAPHICAL DISPLAY OF DATA 
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LOG e NU VS. LOG e RA 







Figure ti.1 lcg Q Nu vs. lcg e Fa (0 degrees roll). 
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LOG e NU 

0.85 1.00 1.05 1.10 1.15 1.20 1.25 



r — 
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LOG e NU VS. LOG e RA 




Fijure B . 2 1 c i] g N vj vs. loj^Ra (10 decrees roll). 
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LOG e NU VS. LOG e RA 




Figure 13.3 1 cg e Nu vs. log^Ra (20 degrees roll). 
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LOG e NU 

0.P5 1.00 1.06 1.10 1.16 1.20 1.25 



LOG e NU VS. LOG e P.A 




I 1 

80 fl.5 7.0 7.5 

log 6 ra 



Fijure n.4 log^Nu vs. lo.jgH.i {JO decrees rolL). 
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LOG e NU 

0.00 1.00 1.00 1.10 1.10 1.20 1.20 



1 



LOG e NU VS. LOG e RA 




Figure B.5 1 o g e M u vs. lcg^Pa (45 degrees roll). 
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LOGpNU 



IX)G e NU = C(0) + LOG e RA 




6.3 8.4 6.6 6.6 6.7 6.8 6.9 7.0 7.1 7 JI 7.6 7.4 



LOG«RA 



O O 

Figure L’.ti Comparison ct J and 2C Results. 
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NUSSELT NO. VS. ROLL ANGLE 




Fi jure C.7 Nu vs. Roll Anqle ( 0 ) . 
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APPENDIX C 



TYPTCAL DATA RUNS 
TABLE C.l 

Typical Results (579 < Ra < 586) 



T = 

s 


296.89 0 = 


0° 








h 


Pi n 


T, ( ° K ) 

D 


T (°K) 
o 


Ra 


Nu 


Nu , 
ch 


(w/m 2 ’°K) 


1 


312.48 


312.25 


579.40 


2.2339 


2.3333 


6 . 162 


2 


312.32 


312.05 


579.31 


2.5850 


2.3283 


7.128 


3 


312 . 50 


312.23 


579.40 


2.9536 


2.3333 


7.043 


M 


312 . 32 


312.03 


572.94 


2.8039 


2.3280 


7.731 


5 


312.37 


312.08 


574.69 


2.7417 


2.3295 


7 .560 


6 


312.59 


312 . 33 


582.43 


2 .4 842 


2.3358 


6.853 


7 


312 .61 


312 . 35 


583.17 


2 .4806 


2 . 3364 


6.843 


8 


312.64 


312 . 34 


583.84 


2.7984 


2.3367 


7.720 


9 


312.41 


312 . 16 


576.73 


2.4066 


2 .3311 


6.637 


10 


312.44 


312 . 16 


577.28 


2 . 7274 


2.3316 


7 .522 


] 1 


312.58 


312.31 


581.97 


2.6458 


2.3354 


7.298 


] 2 


312.66 


312.38 


584.27 


2 .5824 


2.3373 


7 . 124 


] 3 


312. 6 0 


312 . 32 


582.43 


2.6995 


2.3358 


7.446 


1 <i 


312.37 


312.14 


575.71 


2 .2505 


2.3303 


6 .206 


] 5 


312.54 


312.25 


580.31 


2.7625 


2.3340 


7.620 


1 f. 


312.48 


312.19 


578.48 


2.7727 


2.3326 


7.647 


T = 


296.93 0 = 


1D° 










1 


312.76 


312 .47 


586.28 


2.7304 


2 .3389 


7.534 


O 


312.59 


312.33 


581.14 


2.5429 


2.3347 


7.015 


0 


312.78 


312.51 


587.20 


2.5123 


2.3396 


6.932 


1» 


312.59 


312.31 


580.77 


2.7581 


2 . 3 344 


7 .609 


( 

vJ 


312.64 


312.36 


582.51 


2.6972 


2 .3358 


7.441 


6 


312.81 


312.55 


588.37 


2.4534 


2 . 3406 


6.770 


7 


312 .89 


312.63 


590.95 


2.4434 


2 . 3426 


6.743 


8 


312.86 


312.57 


589 .49 


2.7638 


2 . 3415 


7.627 


9 


312.63 


312 .38 


582 .70 


2 . 376 3 


2.3360 


6 .556 


10 


312.67 


312.38 


583.25 


2.6932 


2 .3364 


7.430 


1 ]. 


312.86 


312.58 


589.76 


2.6057 


2.3417 


7 .191 


12 


312.87 


312 . 60 


588.37 


2.5903 


2.3420 


7.038 


1 3 


312 . 82 


312 .54 


590.21 


2 .665 8 


2.3406 


7 . 356 


14 


312.59 


312.31 


580.77 


2 . 7581 


2.3344 


7.609 


1 5 


312 . 76 


312.47 


580.25 


2 .7304 


2.3389 


7 .534 


] b 


312.71 


312.42 


584.44 


2.7380 


2.3374 


7 .554 
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TABLE C.l 



Typical Results (579 < Ra < 586) (cont'd.) 



1' = 296.93 0 = 20° 

s 



] 


312.71 


312.42 


584.44 


2.7380 


2.3374 


7.554 


2 


312.54 


312.27 


579 . 29 


2.5523 


2.3332 


7.040 


i 


312.72 


312.46 


585.36 


2 . 5216 


2.3381 


6.958 




312 . 54 


312 .25 


578.92 


2.7684 


2.3329 


7.636 


h 


312.59 


312.31 


580.68 


2.7047 


2.3343 


7.461 


[. 


312.76 


312.44 


585.84 


2.9965 


2.3384 


8.268 


7 


312.78 


312 . 52 


587.28 


2.4811 


2.3397 


6.791 


8 


312.81 


312.51 


587.65 


2.7738 


2.3400 


7.654 


f j 


312.58 


312.33 


580.86 


2.3851 


2 . 3345 


6.579 


10 


312.56 


312.27 


579.57 


2 . 7130 


2.3343 


7.484 


] 1 


312.75 


312.47 


586.09 


2.6246 


2 . 3387 


7.242 


1 2 


312.82 


312.54 


588.37 


2.5597 


2 . 3406 


7.063 


1 8 


312.77 


312.48 


586.56 


2.6733 


2.3391 


7.376 


1 <i 


312.54 


312.25 


578.92 


2.7680 


2.3329 


7.636 


] 5 


312.59 


312.31 


580.77 


2.7584 


2.3344 


7.609 


1 f> 


312 . 59 


312 . 31 


580.77 


2.7584 


2 .3344 


7.609 


T - 

s 


296.94 0 


= 30° 










1 


312 .65 


312 . 36 


582.81 


2.7480 


2.3359 


7 .581 


O 


312.54 


312.27 


579.29 


2.5523 


2.3332 


7.040 


3 


312.67 


312.40 


583.52 


2.5332 


2.3366 


6.989 




312.54 


312.25 


578.92 


2.7684 


2.3329 


7.636 


i 

K.) 


312.53 


312 .31 


579.76 


2.1768 


2.3336 


6.005 




912.76 


312.44 


585.64 


2.9965 


2.3384 


8.268 


7 


312.78 


312 .52 


587.28 


2.4611 


2 . 3 39 7 


6 . 791 


8 


312.58 


312.27 


579.94 


2.9250 


2.3337 


8.068 


9 


312.56 


312.27 


579.97 


2.7130 


2.3334 


7.484 


10 


312.70 


312.51 


586.73 


2.2509 


2.3392 


6.211 


11 


312.69 


3J2 . 42 


584.27 


2.6842 


2.3373 


7.268 


12 


312.52 


312.24 


588.37 


2.6599 


2 .3406 


7.063 


1 3 


312.77 


312.48 


586.56 


2.6733 


2.3391 


7.376 


1M 


312.48 


312 . 14 


575.25 


2 .7890 


2.3299 


7.692 


1 5 


312 . 59 


312.31 


580.77 


2.7581 


2.3344 


7 . 609 


16 


312.59 


312.31 


580.77 


2.7581 


2.3344 


7.609 



MO 



T 

s 

1 

2 

3 

4 

5 

C 

7 

8 

5 

1 0 

1 1 

1 2 

J 3 

14 

1 8 

1G 



TABLE C.l 



Typical Results (579 < Ra < 586) (cont’d.) 



296.96 0 = 45° 



312 , 


.65 


312 , 


.42 


312 . 


.48 


312 . 


.22 


312 . 


.67 


312 . 


.40 


312 . 


.48 


312 . 


. 19 


312. 


.53 


312 . 


.25 


312 . 


. 76 


312 . 


.49 


312 . 


.78 


312 . 


. 52 


312 . 


. 7 5 


312 . 


.51 


312 . 


.58 


312 . 


.33 


312 . 


.50 


312 . 


.22 


312. 


, G9 


312 . 


. 4 2 


312. 


. 82 


312 . 


.54 


3 12, 


. 71 


312 . 


.43 


312 . 


.43 


312 . 


.14 


312 . 


.59 


312. 


.31 


312 . 


. 59 


312 . 


. 31 



582.12 2.21640 
576.06 2.56722 
582.12 2.53847 
575.69 2.78452 
577.44 2.72293 
585.15 2.46985 
585.89 2.46626 
585.35 2.25553 
679.47 2.39009 
576.34 2.72653 
582.86 2.63975 
586.99 2.56502 
583.31 2.69100 
573.85 2.79496 
579.38 2.76402 
579.38 2.76402 



2 . 33550 


6.115 


2 . 33057 


7.081 


2.33550 


7 .004 


2 . 33024 


7.681 


2.33169 


7.511 


2 .33796 


6 .815 


2.33856 


6.806 


2.33812 


6.224 


2.33335 


6 . 594 


2.33079 


7.521 


2 . 33611 


7.283 


2.33945 


7.078 


2 . 3364 7 


7.425 


2 .32875 


7 .709 


2.33328 


7.625 


2 .33328 


7.625 
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TABLE C. 2 



Typical Results (1305 < Ra < 1380) 



T =317 
s 


.83 0 = 


0° 










Pin 


T b (°K) 


0 

o 

* 
V ✓ 


Ra 


Nu 


Nu ch 


n 

( w/m^ * ' 


1 


338.49 


338.04 


1304 . 18 


2 .97416 


2 . 76085 


6.553 


2 


338.49 


337.99 


1302.95 


3.27379 


2 . 76029 


7 .213 


3 


338 . 52 


338.06 


1305.56 


3.07919 


2.76146 


6.784 


4 


338.44 


337.04 


1300.20 


3.28048 


2 .75905 


7.227 


5 


338.53 


338.05 


1305.56 


3.15390 


2 . 76146 


6.949 


0 


338.49 


338.02 


1303.57 


3.12388 


2.76057 


6 .882 


7 


338.52 


338.05 


1305.40 


3.11699 


2 .76139 


6.868 


8 


338.50 


338.01 


1303.72 


3.23238 


2 .76064 


7.122 


9 


338.53 


338.08 


1306.47 


3.00352 


2 . 76188 


6.618 


10 


338.53 


338.08 


1306.31 


2.96674 


2 .76180 


6.537 


1 1 


338.07 


338.21 


1313.81 


2 .9 8347 


2 .76517 


6 .575 


12 


338.63 


338.17 


1311.66 


3 .06200 


2.76421 


6.747 


13 


3 3 8.50 


333.01 


1303.73 


3.23238 


2 .76064 


7.122 


14 


338.51 


338.04 


1 304 .78 


3.04564 


2 . 76111 


6 .710 


1 5 


338.55 


338.07 


1 306 .6 3 


3.18820 


2 .76195 


7.025 


] G 


338.51 


338.01 


1304 .02 


3 . 30642 


2.76077 


7.285 


T =314 
s 


.72 0 = 


10° 










1 


334.82 


334.35 


1326.67 


3.23972 


2 .77102 


7.075 


0 

L- 


334.87 


3 34 . 42 


1330.20 


3 .076 32 


2.77259 


6.718 


3 


334 .81 


334 .33 


1325.54 


3 .28189 


2 .77052 


7 . 166 


4 


334.79 


334.35 


1325 . 88 


3.04951 


2 .77067 


6.659 


r 

o 


334 .76 


334.32 


1323.93 


3.97872 


2 . 76980 


6.504 


G 


334.82 


334.38 


1327.63 


3 .00764 


2.77145 


6.568 


7 


334.75 


334.28 


1322.49 


3.24195 


2 .76916 


7.101 


8 


334.88 


334 . 39 


1329.72 


3.34543 


2.77238 


7 . 306 


9 


334 . 87 


3 34 . 4 2 


1330.06 


3.03809 


2.77253 


6.635 


10 


334.81 


334.37 


1326 . 84 


3.05689 


2 .77110 


6.654 


11 


334 . 88 


334.41 


1329 . 89 


3.23042 


2 .7724 5 


7.055 


1 2 


334.92 


334.41 


1331 .16 


3.45717 


2.77302 


7 .550 


1 3 


334.95 


334.49 


1 3 34 .06 


3 .06572 


2.77430 


6.696 


1 4 


334.80 


334.34 


1325.88 


3 .12721 


2 .77067 


6 .829 


1 5 


3 34 . 9 3 


334.47 


1333.24 


3 . 18207 


2 .77394 


6 .950 


] 0 


334. 80 


334 .33 


13 2 5 . 5 4 


3.20581 


2.77052 


7 .000 



TABLE C . 2 



Typical Results 



T =311 
s 


.24 


0 = 


10° 


1 


330.59 


330.16 


2 


3 30 . 4 8 


330.07 


3 


330.48 


330 .04 


4 


330.43 


329.98 


r 

o 


330 . 59 


330 . 14 


b 


330.54 


330.07 


7 


330.43 


330 . 00 


8 


330.59 


330.16 


9 


330.54 • 


930.09 


3 0 


330.43 


329.98 


n 


380.69 


330 . 14 


] 2 


330.60 


330.14 


13 


330 . 53 


330.08 


14 


3 30 . 4 8 


330 . 05 


15 


330 . 59 


330 . 16 


1C 


330.37 


329.92 


T =310. 


65 


0 = 


20° 



s 



1 


330 . 


, 32 


329 . 


. 88 


9 


330 , 


. 32 


329 , 


,89 


3 


3 30 , 


.26 


329 , 


. 82 


4 


330 . 


.26 


329 . 


. 82 


9 


3 30 , 


. 32 


329 , 


.87 


6 


330 , 


.26 


329 . 


. 80 


/ 


330 , 


.26 


3 29 . 


.78 


8 


3 30 . 


. 37 


329 . 


. 88 


9 


330 , 


. 4 3 


329 . 


.93 


1 0 


3 30 , 


. 32 


329 . 


.87 


1 1 


330 , 


.43 


329 , 


.98 


1 2 


3 30 , 


.43 


329 , 


.98 


1 3 


330 , 


.42 


T29 . 


. 9 3 


14 


330 . 


. 32 


329 , 


. 82 


15 


330 


.4 3 


329 . 


.9 8 


1 6 


3 30 . 


.26 


329 , 


. 82 



< Ra < 1380) (Cont.) 



11 


3.11643 


2 .78019 


6.736 


81 


3 .01289 


2 .77741 


6.512 


97 


3 .21805 


2 .77704 


6.955 


56 


3.22605 


2.77553 


6.972 


60 


3.23823 


2.77997 


7.000 


69 


3 . 370 37 


2.77824 


7.285 


07 


3.10318 


2.77576 


6 . 706 


11 


3.11653 


2 . 78019 


6.736 


38 


3 . 20669 


2 . 77855 


6.931 


4 6 


3.22605 


2.77553 


6.972 


60 


3.23823 


2 . 77997 


7 .000 


60 


3 . 31969 


2.77997 


7.176 


87 


3 .25068 


2.77832 


7.026 


30 


3.13528 


2 .77719 


6.776 


93 


3 . 15756 


2 .78012 


6 .825 


00 


3.27927 


2 .77396 


7 .086 



21 


3.06971 


2 . 79 335 


6.628 


54 


2 .9 8876 


2 . 79 350 


6.453 


47 


3.16018 


2 .79173 


6.823 


47 


3 . 16018 


2 .79173 


6.823 


70 


3.18968 


2.79313 


6.887 


95 


3.28058 


2 .79151 


7.082 


26 


3.44194 


2 .79121 


7.431 


91 


3.46225 


2 .79409 


7.476 


96 


3.53134 


2 . 79 541 


7.625 


88 


3.14908 


2 .79 321 


6 . 799 


67 


3 .13060 


2 . 79614 


6.760 


85 


3 . 1699 2 


2.69622 


6.845 


82 


3.49360 


2.79535 


7.544 


17 


3.55229 


2 .7924 7 


7.669 


67 


3 .13060 


2.79614 


6 .760 


4 7 


3 .16018 


2.79173 


6.823 



(1305 

1 34 7 . 

1 340 . 

1339 . 

1336 . 

1 346 . 

1 34 2 . 

1337 . 

134 7 . 

1 34 3. 

1 936 . 

1 346 . 

1346 . 

1 342 . 

1 340 . 

134 6 . 

133 3 . 

1377 . 

1 37 7 . 

1 37 3 . 

13 7 3 . 

1 376 . 

137 2 . 

1372 . 

1 378 . 

1381 . 

1 3 7 6 . 

1 38 3 . 

13 8 3 . 

1381 . 

1 37 5 . 

1 38 3 . 

1 373 . 



TABLE C . 2 



Typical Results (1305 



T =303.00 0 = 30° 

5 



1 


319.53 


319.19 


1317 . 


2 


319.53 


319.18 


1318. 


3 


319.56 


319.18 


1319 . 


9 


319.98 


319.09 


1313 . 


5 


319 .55 


319 .18 


1319 . 


6 


319.97 


319.13 


1319 . 


7 


319.51 


319.12 


1315 . 


8 


319.52 


319 . 19 


1316 . 


9 


319 .55 


319.17 


1319 . 


10 


319.99 


319.12 


1315 . 


1 1 


319 .52 


319.17 


1317 . 


12 


319.99 


319.09 


1312 . 


13 


319.56 


319.22 


1321 . 


19 


319.59 


391.17 


1318 . 


15 


319.97 


319.08 


1312 . 


19 


319.57 


319.22 


1321 . 


T =299 
s 


.02 0 = 


30° 




1 


319.53 


319 .19 


1 320 . 


l 


319.53 


319.20 


1 320 . 


3 


3 19.5 7 


319.23 


1 32 3 . 


9 


319.57 


319.29 


1323 . 


1 


3 19.5 6 


319.23 


1 322. 


\ t 


319.97 


319.19 


1316 . 


7 


319.52 


319 .18 


1319 . 


B 


319.56 


319.23 


1 32 3 . 


' J 


319.61 


319.29 


1 32 7 . 


1 n 


319.50 


319.18 


1319 . 


1 1 


319.56 


319.23 


1 32 3 . 


12 


319.51 


319 .19 


1317 . 


1 3 


319.99 


319 . 12 


1316 . 


19 


319 . 9 1 


319 .15 


1318 . 


1 5 


319.97 


319 .10 


1319. 


16 


319.53 


319 .20 


1 320 . 



Ra < 1380) (Cont.) 



3.31536 


2.76727 


6.983 


3.07)08 


2.76789 


6.968 


3.25899 


2 . 76891 


6 .869 


3.37931 


2.76576 


7.106 


3 . 16997 


2.76829 


6.665 


2.98398 


2 .76611 


6.285 


3 . 31993 


2.76657 


6.992 


3.26802 


2.76700 


6.883 


3.26120 


2 . 76 806 


6.869 


3 .27261 


2.76629 


6.892 


2.97396 


2.76753 


6.263 


3.03882 


2.76515 


6.900 


2.96517 


2.76899 


6 .296 


3 .21905 


2.76779 


6.770 


3.92929 


2.76532 


7 .211 


3.06151 


2.76913 


6.999 



3.10591 


2.76886 


6 .962 


3.10715 


2.76909 


6.969 


3 . 20 306 


2.77017 


6.669 


3 .0997 8 


2 .7 7027 


6.335 


3.09922 


2 .769 89 


6 .399 


3.06911 


2 .76706 


6 . 379 


3.10825 


2.76898 


6.966 


3.09906 


2 . 76999 


6.998 


2.93205 


2 .77196 


6 .101 


2.9590s 


2 .76819 


6.196 


3.09589 


2 .77009 


6.337 


3.92918 


2.76752 


7 . 139 


3.98656 


2 . 76686 


7.253 


3 .37939 


2 .767 81 


7 .020 


3.93809 


2 . 766 21 


7.152 


3.10715 


2.76909 


6.969 



< 

29 

68 

86 

92 

98 

72 

73 

70 

06 

10 

88 

56 

09 

97 

95 

9 6 

50 

91 

9 9 

66 

81 

97 

69 

03 

9 6 

00 

2 5 

51 

02 

15 

57 

91 



99 



TABLE C . 2 



Typical Results (1305 < Ra < 1380) (Cont.) 



T =301 . 86 0 

s 



= 30.0° 



1 


318.91 


318.58 


1378.97 


2.79324 


2 . 7944 5 


5 . 869 


2 


318 . 87 


318.48 


1373.41 


3 .27655 


2 . 79204 


6.884 


3 


3)8.96 


318 . 56 


1379.78 


3.35353 


2 . 794 80 


7.047 


4 


318 . 9 6 


318.57 


1379 .99 


3.21071 


2 . 794 89 


6.747 


r 


319.00 


318.62 


1383.17 


3 .2021 3 


2 . 79627 


6.729 


6 


318.97 


318.59 


1381.18 


3 .20759 


2 . 79540 


6 . 740 


7 


319.02 


318.63 


1384.15 


3.24665 


2 . 79669 


6.823 


8 


318.95 


318.56 


1379.36 


3.25983 


2 . 79462 


6.850 


9 


318.89 


318.52 


1375.58 


3 .1296 8 


2.79298 


6 .576 


1 0 


318.94 


318.57 


1379.57 


3.11910 


2 .79471 


6.554 


n 


318.95 


318.57 


1379.57 


3.21197 


2 .79471 


6.749 


] 2 


319 .01 


318.58 


1382 .16 


3.53414 


2.79583 


7.426 


13 


318.99 


318 . 56 


1380 . 76 


3.58582 


2.79522 


7.535 


14 


318.96 


318 . 59 


1380 .76 


3.02125 


2.79522 


6.349 


15 


318.91 


318 . 54 


1377.40 


3.07746 


2 .79377 


6.466 


16 


318.98 


318 . 59 


1381 . 35 


3.25438 


2.79548 


6.838 


T =309 
s 


.50 0 = 


4 5° 










1 


328.21 


327 . 79 


1339.79 


3.07511 


2 .7770 3 


6.6]] 


2 


328.26 


327.83 


1 342 . 74 


3.19372 


2 .778 3 3 


6 . 866 


3 


328.27 


327.87 


1343.98 


2.97942 


2.77887 


6.406 


4 


328.16 


327.72 


1335.93 


3.33687 


2.77532 


7.173 


5 


328.15 


327 . 75 


1336.61 


3.00068 


2 .77562 


6.450 


6 


328.21 


327.84 


1341 . 36 


2.69455 


2 .77772 


5.793 


7 


328.31 


327.84 


134 1 . 36 


2 .77899 


2 .77772 


5.974 


8 


328.31 


327.93 


1347 .4 8 


2 .804 15 


2 .7804 2 


6 .029 


9 


328.21 


327.71 


1337.14 


3 . 71214 


2.77585 


7.980 


10 


328.09 


327 .71 


1333.65 


2.88105 


2.77431 


6.193 


1 1 


328.26 


327.79 


1341 . 36 


3.44784 


2.77772 


7.412 


1 2 


328.26 


327.79 


1 341.54 


3 .4 8963 


2 .77780 


7.502 


1 3 


328.19 


327.69 


1336 .27 


3.75752 


2 .77547 


8.077 


14 


328.11 


327.72 


1334 .52 


2.92132 


2 .77469 


6.279 


15 


328.26 


327.83 


1342.59 


3.23456 


2.77826 


6.954 


1 6 


328.21 


327.75 


1338.55 


3.45635 


2 . 77648 


7.4 30 
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APPENDIX D 



DATA REDUCTION PROGRAM 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



x * * * * * x * -x- -x -x- x- x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x 

* PROGRAM TO CALCULATE HEAT TRANSFER COEFFICIENT ON A * 

* PINNED CIRCULAR DISK * 

X X X "X" 'X- 'X- 'X -X* X "X" X 'X* -X -X- -X* "X' X X" X *X 'X' X *X -X' X -X -X* -X X X X X X X X X X' X XXX X X X' X X X X xxxxxxxxxx 



A 

B 

c 

CB 

CO 

D 

E 

G 

GR 

H 

KA 

KP 

M 

NU 

NUC 

NUAVG 

NUT 

PHI 

PR 

R 

RA 

RAT 

RAV/ 

TB 

TO 

TBF 

TOF 

TB ( I ) 

TO ( I ) 



* -x- x- x x x -x- *DEF I N I T I ON OF VARIABLES * * * * * * * * * * * * * * * * * * 
FACTOR IN CHURCHILL AND CHU CORRELATION USED HERE 
TO SIMPLIFY AND SHORTEN THE EXPRESSION 
LENGTH OF THE SPINE 

SIMPLIFYING FACTOR FOR CHURCHILL AND CHU CORRELA- 
TION 

CORRECTION TO THE SPINE BASE TEMPERATURE 
CORRECTION TO THE SPINE TIP TEMPERATURE 
DIAMETER OF SPINE 

SIMPLIFYING FACTOR FOR THE CHURCHILL AND CHU COR- 
RELATION 

GRAVITATIONAL ACCELERATION 
GRASHOF NUMBER 

HEAT TRANSFER COEFFICIENT (EXPERIMENTAL) 

THERMAL CONDUCTIVITY OF THE FLUID (AIR) 

THERMAL CONDUCTIVITY OF THE SPINE MATERIAL 

A HEAT TRANSFER PROPERTY 

NUSSELT NUMBER (FROM EXPERIMENTAL DATA) 

NUSSELT NUMBER (FROM CHURCHILL AND CHU CORRELA- 
TION) 

AVERAGE NUSSELT NUMBER (FROM EXPERIMENTAL DATA) 
TOTAL OF NUSSELT NUMBERS FROM 16 SPINES (1RUN) 
ROLL ANGLE (DEGREES) 

PRANDTL NUMBER (EXPERIMENTAL) 

RATIO OF TEMPERATURE EXCESSES ( = THB/THO) 
RAYLEIGH NUMBER (EXPERIMENTAL) 

TOTAL RAYLEIGH NUMBER FOR 16 SPINES (1 RUN) 
AVERAGE RAYLEIGH NUMBER (EXPERIMENTAL) 

BASE TEMPERATURE (DEGREES F) UNCORRECTED 
TIP TEMPERATURE (DEGREES F) UNCORRECTED 
CORRECTED BASE TEMP. (DEGREES F) 

CORRECTED TIP TEMP. (DEGREES F) 

CORRECTED BASE TEMPERATURE (KELVIN) 

CORRECTED TIP TEMPERATURE (KELVIN) 



C THB = TEMPERATURE EXCESS AT THE BASE OF SPINE 

C THO = TEMPERATURE EXCESS AT THE TIP OF SPINE 
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oooooo ooo ooo ooo ooooo o o oooo 



TF 

TS 



TEMPERATURE OF FILM BETWEEN SPINE AND FLUID 
TEMPERATURE OF SURROUNDING FLUID (AIR) 



* * * * * * * * **#***#*DECLARATION of VARIABLES * * * * * * * * * * * * * * * * 
REAL D , KP , B , TS , R , M , H , PHI , NU , THB , THO , RAT , RAV , NUT , NUAVG , 
REAL A,C,E,G,TF,VA, ALFA , KA , GR , PR , RA , NUC 
REAL TB(25),T0(25),CB(25),C0(25) 

INTEGER I 



* * * * * * * * * * * * '"'“'INITIALIZATION OF VARIABLES* 
D = 0.0127 
B = 0.0508 
KP = 120 
G = 9.807 
TS = 87.53 
PHI = 45 
NUT = 0 
RAT = 0 



y. y. j 



-X- -X 



■>'. -V- « 



CALCULATE TS IN KELVIN 

TS = (TS - 32) *5/9 + 273.1 5 

PRINT HEADINGS 
WRITE(6,100) TS.PHI 
WRITE(6,1 50) 

READ IN AND CALCULATE CORRECTED TEMPERATURES, DEGREES F 
DO 50 I = 1,16 

READ(5,200)TB(I) , TO ( I ) ,CB(I) ,00(1) 

TBF = TB ( I ) + CB ( I ) 

TOF = TO ( I ) + C0(I) 

CONVERT CORRECTED BASE & TIP TEMP. TO KELVIN 

TB ( I ) = (TBF - 32 ) *5/9 +273.15 
TO ( I ) = (TOF - 32 ) *5/9 +273.15 

CALCULATE THETA BASE AND THETA ZERO (TEMP. EXCESSES) 

THB = TB ( I ) - TS 
THO = TO ( I ) - TS 

CALCULATE RATIO OF THETA B TO THETA 0 (=cosh(MB)) 

R = THB/ THO 
CALCULATE M 

M = ( ALOG (R + SQRT ( (R**2 ) -1 ) ) ) /B 
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C CALCULATE HEAT TRANSFER COEFFICIENT H 
C 

H = (M**2)*KP#D/4 
C 

C CALCULATE FILM TEMPERATURE 
C 

TF = (1 . /4. )*(TB(I) + TO ( I ) ) + (1./2.)*TS 
C 

C CALCULATE PROPERTIES OF AIR AT TF 
C 

KA = 0.004033 + (7.30E-5)*TF 

VA = (-1.22567 + (9.350E-3)*TF)*(1 .OE-5) 

ALFA = (-1.82333 + (0. 01 350*TF) ) *(1 • OE-5 ) 

C 

C CALCULATE NUSSELT NUMBER BASED ON EXPERIMENTAL H 
C 

NU = H*D/KA 
NUT = NUT + NU 
C 

C CALCULATE GRASHOF, PRANDTL , RAYLEIGH, AND NUSSELT NUM- 
C BERS BASED ON CHURCHILL & CHU CORRELATION 
C 

PR = VA/ALFA 

GR = (G*(1 .0/TF)*(TF-TS)*D**3)/VA**2 
RA = GR * PR 
A = 0.387*RA**(1 .0/6.0) 

C = (1.0 + (0. 559/PR)**(9. 0/16. 0))**(8. 0/27.0) 

E = A/C 

NUC = (0.60 + E)**2 
RAT = RAT + RA 
C PRINT OUTPUT 
C 

WRITE(6,250)I,TB(I) , TO ( I ) ,RA,NU,NUC,H 
50 CONTINUE 
C 

RAV = RAT/16.0 
NUAVG = NUT/16.0 
WRITE ( 6 , 300 ) RAV , NUAVG 
C 

STOP 

1 00 FORMAT ( ' 1 ' ,1 5X, 'TS= ' ,F6.2,3X, 'PHI= ' ,F4.1 ) 

1 50 FORMAT (1 7X , 'PIN' , 4X , ' TB ' , 7X , ' TO ' , 7X , ' RA ' , 8X , ' NU ' , 8X , 
ff 'NUC' ,7X, ' H ' ) 

200 FORMAT ( 4 (F 7.3) ) 

250 FORMAT (1 8X, 12 , 2X , F6. 2 , 3X , F6. 2 , 2X , F7 . 2 , 2X , F8 . 5 , 2X, F8 . 5 , 
§ 3X.F6.3) 

300 FORMAT (1 6X, 'RAV = ' , F7 . 2 , 2X , ' NUAVG =',F6.3) 

C 

END 
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